DESY 99-022 
MZ-TH/99-03 



Supersymmetry and CP Violating 
Asymmetries in B^ s Decays 

Gabriela Barenboim a,x and Martti Raidal b ' 2 

& Institut fur Physik, Johannes Gutenberg Universitat, D-55099, Mainz, Germany 
b Theory Group, DESY, Notkestrafle 85, D-22603 Hamburg, Germany 

Abstract 

We study possible effects of supersymmetry (SUSY) in CP asymmetries in non- 
leptonic Bd )S decays in a variety of SUSY flavour models considered in literature. 
We use both mass insertion and vertex mixing methods to calculate squark-gluino 
box diagrams contribution to Bd, s -Bd,s mixings. With the squark mixing parameter 
n = 0.22, and with large new CP phases, it turns out that the CP asymmetries 
to be measured in upcoming B-factories, HERA-B and LHC-B, can be completely 
dominated by the SUSY contribution in almost every considered model. Discrimi- 
nation between the different models can be done by comparing experimental results 
in different decay modes. In some models squark masses up to ~ 5 TeV can be 
probed through these experiments provided the SUSY contribution to B-B mixing 
is at 10% level, \Mf^ V /Mf 2 > | ~ 0.1. This implies that models with heavy squarks 
have a fair chance to be tested in the future CP experiments before LHC. 



1 Introduction 



There are two major questions to be answered in particle physics. One is the possible origin 
of CP violation, to be tested in forthcoming B-factories and the dedicated experiments 
LHC-B and HERA-B. The other one is the possible existence of supersymmetry (SUSY) 
as evidenced by the continuing efforts of both experimentalists and theorists in searching 
for new physics beyond the Standard Model (SM). The SM has specific predictions [1] on 
the size as well as on the patterns of CP violation in Bd, s meson decays which, if disproved 
in future experiments, would signal doubtless the existence of new physics [2]. 
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SUSY theories have been subject of an extensive study from the flavour physics point of 
view since they naturally offer new sources of flavour changing neutral currents (FCNC) 
and CP violation, arising from the Yukawa and SUSY breaking sectors of these models. 
These extra contributions are already constrained in a rather stringent way by the ex- 
perimental bounds on ex, K, B, D mixings, electric dipole moments of electron d e and 
neutron d n as well as on the branching ratio of b — > S7. Many works on radiatively induced 
FCNC processes in SUSY [3] have concentrated on the minimal supersymmetric extension 
of the SM (MSSM), and assumed often the minimal supergravity (mSUGRA) framework 
for the SUSY breaking sector. Updated analyses of mSUGRA imply [4] that AM Bd , ex 
as well as the branching ratio of the processes b — > sll may be enhanced at most by a 
few tens of percent while the branching ratios of b — > svv and K — > -nvv may be reduced 
by only a few percent. The new SUSY CP phases have negligible effects on the decays 
B — > X s 7 and B — > X s ll in mSUGRA. In particular, there is no new large phase shift 
in the B-B mixing [5] implying no major deviations in any CP asymmetry as compared 
with the SM predictions. 

This conclusion changes dramatically as soon as one considers SUSY GUTs or other SUSY 
models without universality in the soft breaking sector. As the squark mass matrices are in 
general independent of the quark mass matrices, they introduce a major source of FCNC 
which, together with new CP phases, can substantially modify the SM predictions for 
CP violating observables to be measured in B-factories and LHC-B. There is a plethora 
of SUSY models, to be considered below, which satisfy all the present phenomenological 
constraints and allow large signals of new physics in these experiments^} 

Many of the recent works on SUSY signatures in CP violating decays of Bd, s mesons 
have concentrated on studying new physics contributions to decay amplitudes. These 
analyses cover the full spectrum of SUSY models. Large deviations from the SM in CP 
asymmetries in b — > S7 are predicted in Ref. [7] and in penguin dominated non-leptonic b 
decays in Ref. [8]. As the new physics contributions to decay amplitudes are non-universal, 
comparison of CP asymmetries in different decay modes allows to find new physics in a 
model independent way. Surprisingly, while the effects of new physics contribution to B- 
B mixing in CP violating lepton asymmetries have been analyzed in almost all possible 
SUSY flavour models in both B^, B s decays [9], the well known mixing effects in CP 
asymmetries in non-leptonic Bd, B s decays have been considered in detail only in a few 
particular models of supersymmetry [10,11]. The mass insertion method has been used to 
describe the internal squark effects in these works. These recent analyses have shown [12] 
that LEP2 experiments pose a serious naturalness problem to all conventional mSUGRA 
models, motivating the search for different SUSY models. In addition, motivated by the 
recent claim [13] that the CP violating observable ex may have a fully supersymmetric 
origin in general SUSY models, we feel encouraged to study the predictions of the full 
range of SUSY flavour models considered in literature on the CP asymmetries to be 
measured in the upcoming experiments. 



3 For a review of SUSY CP violating signals at future colliders see Ref. [6] . 
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Following the model independent analyses of Ref. [14] we study in this letter the CP vio- 
lating asymmetries in B^ s meson decays in the SUSY flavour models classified in Ref. [9] . 
This includes models where FCNC and CP problems are solved using alignment [15-18], 
non-abelian symmetries [18-26] and decoupling of squarks of the first two generations 
[23,27-32]. These cover a large variety of general SUSY flavour models considered in liter- 
ature. We use both the mass insertion method [33,34] as well as the vertex mixing method 
[35] to describe the dominant gluino-squark box diagram's contribution to the meson mix- 
ing. We consider both and B s decays, as they provide complementary information for 
distinguishing new physics from the SM as well as discriminating between different SUSY 
models. We show that despite the quite strong bounds on the squark masses in these 
models, given the predictions for squark mixings, the SUSY CP violating contribution in 
Bd, s systems might dominate over the SM one in almost all the models. In some models 
measurable effects can be achieved for squarks masses as high as several TeV. This implies 
the possibility of discovering SUSY in B-factories before LHC. 



2 New physics in B-B mixing 



Possible ways of evidencing the existence of new physics in B-factory experiments have 
been studied in the literature in a model independent way. Signals of new physics arising 
form new contributions in the decay amplitudes have been studied in [8] and are beyond 
the scope of this paper. If new physics contributes to the B-B mixing, as it is the case with 
the general SUSY models, it has been shown, that the planned experiments will be able to 
distinguish the effects of new CP phases from the effects of the SM Cabibbo-Kobayashi- 
Maskawa (CKM) phase [36]. This can be done by comparing the CP asymmetries in the 
different decay modes of both B d and B s mesons. This important result implies that 
future experiments will allow us to discriminate between different SUSY models, as will 
be shown below. 

In terms of the off-diagonal elements of the 2x2 B^-B^, q — d, s, mixing Hamiltonian 
M — iT/2 the B®-B® mixing phase 4> 9 M is expressed as 



M 12 ~ 2 1 12 m 

\ m? 2 - in 2 • 1 > 

In the analysis that follows, we make use of the fact that M 12 can be parametrized as [14] 

M 12 = Mf 2 M + M™ SY = Mf 2 M (l + he w ) , (2) 

where M^ SY denotes the new SUSY contribution and h = \M^ SY /M^\. The relative 
phase 9 = arg(M 1 5 2 c/5y /Mf 2 M ) is in principle not restricted because it results from the 
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possibly large new phases arising in the squark sector. As | T 12 1 <C |Mi 2 | and Ti 2 ~ Tf 2 M , 
the B°-B° mixing gets modified by the phase 0m 



In the Bd system the new physics contribution is constrained by the measured B\-B\ 
mass difference which may by expressed as AM^ = 2|Mf 2 |. 

The phase (3) directly modifies the CP asymmetries as discussed in [36,14]. Its effects are 
universal in all decay modes. Because the SM predicts almost vanishing CP asymmetries 
in most of B s decays, non-zero experimental results in these measurements unambiguously 
measure the new physics phase 0m- However, in B d system (to be experimentally probed 
in B-factories) the CP asymmetries can be also large within the SM, therefore finding new 
physics contributions to B-B mixing requires an accurate knowledge of the SM predictions 
for these asymmetries. In the presence of new phases, the CP asymmetries which measure 
the SM CKM angle (3 (e.g. in B d — > J/ipK) are modified as acp = — sin2(/9 + 0m) while 
the CP asymmetries which measure the SM CKM angle a (e.g. in Bd — > 7r + 7r~) receive a 
new contribution with the opposite sign, acp = — sin 2 (a — 0m)- As the new phase, 0m, 
cancels out in a + (3 measurements of the third CKM angle 7 are crucial. 

However, if the new physics contribution is large enough so that the measured CP asym- 
metries are outside the allowed SM regions then, the new physics can be traced off un- 
ambiguously also in B-factories. Fortunately, recent global analyses [37] show that the 
precision reached in constraining the CKM matrix is quite good even now, yielding the 
result sin 2/3 = 0.73±0.21 at 95% confidence level. Therefore it follows from Eq. (3) that a 
SUSY contribution of h — 0.1, together with large phase 6, implies measurable deviations 
from the SM. In the following we assume that the minimal detectable value of h is 0.1. 
Note that this assumption is somehow conservative, with better experimental precision, 
in the future much smaller effects can be probed. 



3 SUSY flavour models 

There are two distinct sources of flavour violation in general SUSY models in addition 
to the usual CKM mixing in the quark sector: (i) flavour violating interactions of top 
quark with charged Higgs, and (ii) misalignment between the fermion and sfermion mass 
matrices. The former possibility has been studied in [3] and it is quite constrained to 
lead to significant deviations from the SM. In this work we concentrate on the latter 
possibility. In this case the new physics contribution to the Bd tS mixing is dominated by 
the box diagrams with gluinos, g, and squarks, q, running in the loop. The new flavour 
mixings in the 6x6 down squark mass matrix, 




(3) 
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( Mf L Md2 



LR 



M d R \ Mf R 



(4) 



lead to new flavour changing interactions, which, together with the new CP phases, may 
significantly contribute to Bd yS meson mixings. In Eq. (4), the phenomenological con- 
straints on the off-diagonal squark mixings in 3 x 3 matrices Mf R RL are more stringent 
than in M^i rr [34], and we therefore neglect Mf^ RRL in the following. 

The most oftenly used parameterization of squark mixing effects is called the mass inser- 
tion method [33,34]. It assumes a flavour diagonal gqq vertex and also flavour diagonal 
quark mass matrices and places all the mixing effects, described by the dimensionless 
parameter (S^mn, M,N = L, R, (see Eq. (2.35) in Ref. [9]), in the squark propagators. 
This method is the appropriate one when the squark masses are nearly degenerate but the 
price to be paid for using it, is the introduction of an average squark mass which should be 
chosen in an appropriate way [34]. The second method, the vertex mixing method, deals 
with the mass eigenstates of quarks and squarks with off-diagonal gqq couplings, and con- 
siders only the contribution from the lightest squark generation. Thus the mixing effects 
in this formalism can be characterized by the matrices Kf = VfV^ and K R = V R V R , 
where V£ R and V[ R are the matrices diagonalizing the quark and squark mass matrices, 
respectively. Obviously, this is a better approximation when one generation is much lighter 
than the others. In our numerical calculations we will use both methods depending upon 
which one is the most appropriate for the concrete application to the case at stake. 

The SUSY flavour models were already classified and studied in a comprehensive work, [9] . 
There are three basic mechanisms to suppress the dangerous FCNC. The first one is the 
alignment of squark mass matrices along to the quark ones, so that the mixing matrices 
Kf R are close to unity [15-18]. Motivated by the different behavior of the third family as 
compared to the first two ones, there have been proposed models with non-abelian flavour 
symmetries [18-26]. The approximate flavour symmetries are broken by a small factor 
r) which in our numerical estimates is taken to be equal to the Wolfenstein parameter, 
i] = 0.22. Finally, there are models with super-heavy squarks in the first two generations 
[23,27-32]. To suppress FCNC completely, several of the methods described above may 
be combined in a single realistic model. It has been pointed out that in order to avoid 
fine tuning in the electroweak symmetry breaking [27], and to ensure the positivity of the 
stop mass matrix squared at weak scale [38], there should be bounds on the squark mass 
parameters in consistent heavy squarks models. These, in turn, constrain the most popular 
scenarios of generating the multi TeV masses for the first two generation squarks using 
horizontal U(l) symmetries [28,30,31]. On the other hand, recently it has been proposed 
that all the SUSY soft masses might be at multi TeV scale and the lightness of the third 
generation squarks should be then generated radiatively [39]. In this scenario the right 
bottom squark remains heavy with mass of several TeV. The analysis of such a model is 
beyond the scope of the present work. We shall take a purely phenomenological approach 
and allow squark masses, in particular sbottom masses, to vary from (9(100) GeV to 
several TeV and, study the sensitivity of the CP asymmetries to the squark masses. 
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Table 1 



SUSY models which solve the FCNC problem with A-alignment, B-non-abelian symmetries 
and C-heavy squarks in first two generations. The mixing parameter 77 in LL and RR squark 
mass matrices is taken to be rj = 0.22 in our numerical estimates. The lower bounds on nig 
coming from the measurement of AMs d are presented for each model. The maximal values of 
h = \M?y2 SY /My2 i \ in and B s systems for each model are calculated for the given lower 
bound on rriq. 

It is important to notice that from all the models listed above, some either are not able 
to pass all the phenomenological constraints set by ex, AMx, d e and d n or do not have 
specific prediction for the squark mass matrix textures (for discussion see Ref. [9]). The 
models which have interesting predictions for B-physics, are summarized in Table 1. In 
the first column of this table, we specify the appropriate method for the estimation of 
the SUSY box contribution in each model, while in the second column, we specify the 
mechanism used to suppress FCNC. For a detailed description of each of these models, 
we refer the reader to the original literature which is listed in the third column. The 
predictions for the (23) and (13) flavour mixings in the mass matrices Mi L:RR are also 
given for each model. 



4 Numerical results 



In the SM the quantity Mf 2 M has been calculated including NLO QCD corrections (for 
references see [1]) and reads 

M S 12 M = ^VQC D B B J 2 M Bq M 2 w (V tq V t t) 2 S (z t ), (5) 
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where 



Azt - llz 2 t + z\ 3z f 3 ln* t 

5o( ^ " 4(i - Zt y ~ 2(T^F ' (6) 



with z t = m 2 /m 2 v . Recent updated values for the bag parameter, B Bd , and the decay 
constant, f Bd , are B Bd = 1.29 ±0.08 ±0.06 and f Bd = 175 ±25 MeV, respectively, and the 
QCD correction factor t]qcd takes the value i]qcd — 0.55 ± 0.01 [1]. The meson masses 
we have used for the numerical estimation are, M Bd = 5.279 GeV, M Bs = 5.369 GeV, 
and the value of the top quark mass, m t (m t ) = 165 GeV, is taken from [1]. For the CKM 
mixing elements we use the following values |V^V^| = 0.0084 and |Vt s ^| = 0.040. 

The dominant contribution to Mf^ SY comes from the AB = 2 box diagrams with q, 
g running in the loop. This can be calculated either in the scenario of vertex mixing 
[35] or using mass insertion [33,34]. The latter method is widely used in literature while 
the former, which should be more suitable for models with super heavy squarks, is not. 
The LO [40] and NLO [41] QCD corrections to AF = 2 processes in SUSY models are 
calculated using mass insertion method. There are no calculations of QCD corrections 
using the vertex mixing method. Since one of the aims of this paper is to compare the 
results in vertex mixing and mass insertion we neglect QCD corrections here and work at 
the level of electroweak box diagrams. Because the QCD corrections are known to enhance 
the SUSY contribution significantly (for the Wilson coefficients using mass insertion see, 
e.g., Ref. [11]), the results in this work should be regarded as a conservative estimate of 
the SUSY contribution to CP asymmetries. Moreover, as what matters in CP asymmetries 
are the ratios of amplitudes, QCD corrections tend to cancel. 

Within the vertex mixing approximation, M^ SY is given by [35] 
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where x — iri~/m^ : % — 1, 2 for B d and B s , respectively and we take (K L R ) 33 ~ 1. Within 
the mass insertion notation, Mf^ SY takes the form [34] 
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For our numerical estimates we use a s (m b ) = 0.222, rrib(mb) = 4.4 GeV and x = m 2 ,/m'~ = 
1, unless stated otherwise. 

To begin with, we have to take into account the constraint on the squark masses coming 
from the measurement AM# d = 0.470 ± 0.019 ps" 1 . In order to do so we require that 
our calculated AM Bd = 2|iWf 2 | coincides, within errors, with the experimental value. The 
errors are completely dominated by the errors of Bs d and fs d . The lower bounds on the 
squark masses, m 5 -, for rj = 0.22 and for a fixed value of 8, 9 = n/2, (that takes into 
account the case where the SUSY CP violation is large) for each model are presented 
in Table 1. The interpretation of each of these numbers depends on the method used 
to calculate the SUSY contribution. For the mass insertion method (squark masses are 
nearly degenerate) this is the average squark mass, while for the vertex mixing method 
(squarks of the first two generation are decoupled) it bounds the sbottom masses. 

In principle, one can get bounds on the squark masses and mixings also from the measure- 
ment of b — > 57 branching ratio. However, as we neglect the LR mixings in Eq. (4) here, 
the bound on the (23) squark mixing in LL is at best of order 0(1) for rriq = 500 GeV 
[34], even assuming that the SUSY contribution to the branching ratio is less than 10%. 
Comparing this bound with the model predictions in Table 1 we conclude that b — > 57 
does not impose constraints on our results. 

Comparing the squark mass bounds in Table 1, it follows that in models with the same LL 
and RR mixings the vertex mixing method gives much larger contribution to M12 than 
the mass insertion one (compare models [15], [20], [23] b with [23] a). Since the models 
we have considered here have very different predictions for the squark mixings the mass 
bounds vary from the direct Tevatron bound rriq ^ 260 GeV to 1.4 TeV. 

For the given squark mass bounds, the maximum allowed values of h = |M£> V /M^\ 
for both B^ s systems in each model are also given in Table 1. For the models [16] and 
[18] a, as there is no bound coming from AM Bd and in order to be conservative, we give 
h in the B s system for = 1.4 TeV. As it is evident from the results shown in the table, 
the SUSY contribution dominates in all the models besides model [17], where it can be 
still measurable in the Bd system. As each model predicts different squark mixings the 
maximum h is different in the B^ and B s systems. 

The largest SUSY contribution to the B d -B d mixing is predicted within models [18] b and 
[23] a, and that to the B s -B s mixing by the models [16], [18] a and b, and [23] a. To study 
how large squark masses can possibly be probed in the prospective CP experiments we 
plot in Fig. 1 the dependence of h on rriq in these models (Bd in Fig. 1 (a) and B s Fig. 
1 (b)). The solid curve in Fig. 1 (a) is for model [18] b while in Fig. 1 (b) it is for all 
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(a) (b) 

Fig. 1. h = \Mf 2 USY /Mf 2 M \ as a function of m g - (in GeV) for x = 1 for B d (figure (a)) and B s 
(figure (b)) systems. The solid curves are for the model [18] b (figure (a)) and for the models 
[16], [18] a,b (figure (b)). The dotted curves denote always models [20], [23] b, the long dashed 
curves models [29], [31], and the short dashed curves model [23] a. 




x 



(a) (b) 

Fig. 2. h = \Alf 2 USY /M? 2 M \ as a function of x = rfig/nig for nig = 2.5 TeV for B d (figure (a)) 
and B s (figure (b)) systems. The notation is the same as in Fig. 1. 



the models in [16,18] because their predictions coincide. The dotted curves denote always 
models [20], [23] b, the long dashed curves models [29], [31], and the short dashed curves 
model [23] a. Numerically h is sizable, i.e. it has values of order 0.1, up to squark masses 
of several TeV as can be seen from the graphs. 

To account for the case with very heavy squarks we plot, in Fig. 2, h as a function of x for 
a fixed value rriq = 2.5 TeV. The models and notation are the same as in Fig. 1. Because 
of the large squark mixings, h is always largest in model [18] b. However, in models where 
the first two squark generations are decoupled (models [23] a, [29,31] ) and where the 
vertex mixing method is used, the value of h is in general increased for small values of x. 
Therefore, models with large sbottom masses (gluinos are expected to be light) can still 
be tested in future experiments. 
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5 Discussion and conclusions 



As discussed in Section 2, the measurements of non-leptonic CP asymmetries through 
different Bd and B s decay modes, allow one to discriminate between new physics and 
the SM description of CP violation arising solely from the CKM phase. Our studies 
show that these measurements can also discriminate between the different SUSY models. 
According to Eq. (3), it follows from Table 1, that SUSY can completely dominate the 
CP asymmetries to be measured in B-factories and LHC-B. This is so even once the very 
strong bounds on the squark masses in Table 1 are properly taken into account. Large 
effects in non-leptonic CP asymmetries in Bd system can be expected almost for any 
prediction of the relative size of the LL and RR mixings. This should be compared with 
the leptonic asymmetries where new physics can be probed only if LL = RR [9]. As the 
different models, in general, have different predictions for the (13) and (23) mixings of the 
squark mass matrices (which may differ from the CKM mixing structure), the comparison 
of the CP asymmetries in Bd and B s allows a clear discrimination between models. For 
example, models [16] and [18] a can give large contributions to CP asymmetries only in 
B s decays while the rest of the models can modify either B d or both mesons decays. 

It is also instructive to compare the mass insertion method with the vertex mixing one. It 
follows from Table 1 that in models with super heavy squarks in the first two generations, 
the SUSY contribution for the same value of rrig is larger than in other models. This is 
because in models with nearly degenerate squarks, the GIM mechanism is operative while 
in heavy squark models, only the b squarks contribute. Also, the results in the vertex 
mixing case are much more independent of the relative magnitude of LL and RR mixings 
(compare, for example, models [29,31] with models [21,25]) implying that the new physics 
should affect the CP asymmetries in super heavy squark models more strongly. 

Finally, it remains to be answered how large the squark masses can be and still give 
measurable effects in CP asymmetries. This question becomes very interesting in light of 
the results of Ref. [12]. These authors claim that LEP2 results may indicate that colored 
particles have masses of few TeV. Even in a very disadvantageous scenario, where the 
SUSY effects can be as low as 10%, i.e. h — 0.1, but provided that the phase 9 is large 
enough, there will be still observable effects in upcoming B-factories as shown in Section 
2. As can be seen from Fig. 1 (a), rriq < 5.1 TeV and m~ b < 3 TeV in models [18] b and 
[23] a, respectively, can be probed for 77 = 0.22. The B s decays (Fig. 1 (b)) are somehow 
less sensitive, implying rriq ^ 4.5 TeV and m~ b ^ 2.6 TeV in models [16], [18] a,b and [23] 
a, respectively. If the sensitivity of the future experiments would turn out to be better 
than 10% then, higher masses can possibly be probed. For smaller values of rj the mass 
reach scales linearly. These results have important and far reaching implications on the 
heavy squark models. As discussed, in some models the right bottom squarks might have 
masses of several TeV. Have the new large SUSY CP phases happen in this sector, our 
results would imply (plots in Fig. 2) that TeV masses can still give observable effects in 
B-factories. 
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In conclusion, we have shown that, with the squark mixing parameter 77 = 0.22 and with 
large new SUSY CP phases, the CP asymmetry measurements in upcoming B-factories, 
HERA-B and LHC-B can be completely dominated by the SUSY contribution in almost 
every SUSY flavour model that we have considered. Discrimination between the different 
models can be done by comparing the results of the different decay modes. Assuming that 
SUSY effects at the level of 10% are still measurable, namely a \M?^ SY /M? 2 M \ ~ 0.1 can 
be tested, in some models the sensitivity is enough to explore squark masses up to ~ 5 
TeV. This implies that the models with heavy squarks have a great chance of being tested 
in future experiments. 
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